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Abstract: Land surface temperature (LST) is one of the most important variables measured 
by satellite remote sensing. Public domain data are available from the newly operational 
Landsat-8 Thermal Infrared Sensor (TIRS). This paper presents an adjustment of the split 
window algorithm (SWA) for TIRS that uses atmospheric transmittance and land surface 
emissivity (LSE) as inputs. Various alternatives for estimating these SWA inputs are 
reviewed, and a sensitivity analysis of the SWA to misestimating the input parameters is 
performed. The accuracy of the current development was assessed using simulated 
Modtran data. The root mean square error (RMSE) of the simulated LST was calculated as 
0.93 °C. This SWA development is leading to progress in the determination of LST by 
Landsat-8 TIRS. 
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1. Introduction 

Land surface temperature (LST) is related to surface energy and water balance, at local through 
global scales, with principal significance for a wide variety of applications, such as climate change, 
urban climate, the hydrological cycle, and vegetation monitoring [1^]. LST variations in space and 
time, measured by satellite remote sensing, are used for the estimation of a multitude of geophysical 
variables, such as evapotranspiration, vegetation water stress, soil moisture, and thermal inertia [5-7]. 
With the increasing recognition of the importance of LST, methods for its estimation from space have 
continuously been developed [8]. In recent decades, sensors, such as the Moderate -resolution Imaging 
Spectroradiometer (MODIS) and the Advanced Very High Resolution Radiometer (AVHRR), have 
provided public domain global thermal data twice daily, using two long- wave infrared (LWIR) bands. 
Landsat-5 Thematic Mapper (TM) and Landsat-7 Enhanced Thematic Mapper Plus (ETM-i-) provide 
thermal data using just one long- wave infrared (LWIR) band, with a higher spatial resolution but with 
a 16-day temporal resolution. Since satellite remote sensing provides a repetitive synoptic view in 
short intervals of the Earth's surface, it is a vital tool for monitoring LST. 

Landsat-8 was successfully launched on 11 February 2013 and deployed into orbit with two 
instruments on-board: (1) the Operational Land Imager (OLI) with nine spectral bands in the visual 
(VIS), near infrared (NIR), and the shortwave infrared (SWIR) spectral regions; and (2) the Thermal 
Infrared Sensor (TIRS) with two spectral bands in the LWIR. The relative spectral response of the 
TIRS bands is presented in Figure 1. The two TIRS bands were selected to enable the atmospheric 
correction of the thermal data using a spht window algorithm (SWA) [9,10]. The use of two separate, 
relatively narrow, thermal bands has been shown to minimize the error in the retrieval of LST [11]. 
The spatial resolution of TIRS data is 100 m with a revisit time of 16 days, and as a result, applications 
are different than those of other sensors with coarser spatial resolutions and shorter revisit times. While 
Landsat-8 images are already freely distributed through the U.S. Geological Survey (USGS), to the 
best of our knowledge, no SWA for LST retrieval from TIRS has been published. Although several 
SWAs have been developed for use with other sensors [12-15], some adaptations are required in order 
to implement them for the TIRS spectral bands. Therefore, the objective of this letter is to develop a 
SWA, adapted for use with Landsat-8 TIRS data, along with its accuracy assessment. 

Figure 1. Landsat-8 TIRS bands' relative spectral response functions (the data can be seen in [16]). 
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2. Split Window Algorithm (SWA) 

The SWA was first proposed by McMillin [17] who suggested using the differences in the 
atmospheric absorbance of two adjacent LWIR bands in order to accurately retrieve the sea surface 
temperature (SST). In order to make the transition from SST to LST retrieval, one has to assume the 
land surface emissivity (LSE) in both bands a priori [12]. Qin et al. [18] have presented a SWA for 
AVHRR that requires only two essential variables: LSE and atmospheric transmittance. They tested 
their algorithm and found its accuracy to be 1.75 °C for real world data. Additionally, they found that 
it was preferable to the other SWAs that also performed well but required some parameters that are 
difficult to estimate [19,20]. Therefore, the SWA suggested by Qin et al. [18] was chosen to be 
adapted for TIRS, not only because it was tested and proved to be accurate, but also because it is 
reasonable to estimate its input parameters, as will be discussed in the next section of the paper. 
Further, the current work complements that of Qin et al. [18] to which we refer the reader for 
theoretical background and algorithm development. 

The adaptation of the SWA for TIRS bands relies on the determination of parameter for the 
TIRS-specific spectral bands. L, is defined as: 

Li = Bi(T) /[dBi(T)/dT] (1) 

in which L, has the dimension of temperature in Kelvin, 5,(7) is the Planck function radiance, 
spectrally integrated over each of the TIRS bands at temperature T, and dBi(T) is the derivative of the 
Planck function for band / at temperature T. Therefore, dBi(T)/dT can be calculated as: 

dBi(T)/dT ~ [Bi(T + AT) - Bi(T)]/AT (2) 

We computed L, numerically by using dBi{T)ldT and accurately fit it using a linear regression 
Li = at + biT. For the temperature range of 0-60 °C, Lw = -64.4661 + 0.4398r = 0.9968, standard 
error of estimation (SEE) = 0.1643) and Ln = -68.8678 + 0.4755r (r' =0.9967, SEE = 0.1687) 
(Figure 2). However, coefficients a, and bi vary when computed for different temperature ranges, as 
presented in Table 1. From this table, we can see that as the range of 7 decreases, a better accuracy, or 
a lower SEE, is achieved. In this paper, the SWA accuracy assessment was carried out for the strict 
case of the extreme temperature range of 0-60 °C. However, in order to obtain the LST as accurately 
as possible, it is advisable to select coefficients according to the temperature range in the image. 

Figure 2. Variations of (A) Planck's radiance and (B) L, parameter with temperature for 
each of the TIRS bands. 
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Table 1. Regression coefficients (a, b), coefficient of determination (r ), and Standard 
Error of Estimate (SEE) of the derivative of the Planck function, dB, for Landsat-8 TIRS 
bands 10 and 1 1, for estimating parameter L; at different ranges of temperatures (7). 



T Range (°C) 


aio 


bio 


4. 


SEEio 


an 


bu 




SEEu 


0-30 


-59.1391 


0.4213 


0.9991 


0.0424 


-63.3921 


0.4565 


0.9991 


0.0438 


0-40 


-60.9196 


0.4276 


0.9985 


0.0746 


-65.2240 


0.4629 


0.9985 


0.0769 


10-40 


-62.8065 


0.4338 


0.9992 


0.0415 


-67.1728 


0.4694 


0.9992 


0.0427 


10-50 


-64.6081 


04399 


0.9986 


0.0730 


-69.0215 


0.4756 


0.9986 


0.0750 



The SWA is derived from a first-order Taylor-series linearization of the radiative transfer equation [21], 
and its formulation by Qin et al. [18] takes the general form of: 



(3) 



where Ts is the LST, Tjo and Tuare the brightness temperatures of TIRS bands 10 and 11, respectively, 
and Ao, Aj, and A2 are coefficients determined by the atmospheric transmittance and LSE in both 
TIRS bands: 

Ao = E^a^o + (4a) 
v4i = 1 + i4 + E^b-^Q (4b) 
A2=A + E^bu (4c) 
Based on the algorithm developed by Qin et al. [18], we define: 

Ci = SiTiie) (5) 

Di = [l-Ti(0)][l + (l-£OTi(0)] (6) 

where £, is the LSE of band i and r,(0 is the atmospheric transmittance for a given zenith view angle 6 
in band i . Accordingly, the parameters in Equation (4) are defined: 

A = D^o/Eo (7) 

^1 =£>ii(l-Cio-£'io)/5o (8) 

£2 = Dio(l-Cii-Dii)/£o (9) 

^0 ~ ^11^10 ~ ^10^11 (10) 

The algorithm in Equation (4) uses the estimation of two geophysical parameters, namely 
atmospheric transmittance and LSE, in order to estimate the LST. The next section will briefly discuss 
the estimation of these parameters. 

3. Estimating the Input Parameters of the SWA 

3.1. Determination of Atmospheric Transmittance 

Some SWAs were developed for wide swath sensors (e.g., NOAA-AVHRR, MODIS), and 
consequently emphasize the correction for zenith view angle effects on the atmospheric transmittance. 
However, in the case of the Landsat-8 TIRS at an altitude of 705 km with a swath of 185 km, the 
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maximum zenith view angle is about 7.5°. At that angle, the effect on the atmospheric transmittance in 
both LWIR bands is negligible [18]. Thus, the term 0 may be removed from Equations (5) and (6) for 
the purpose of implementing this SWA for TIRS. 

In addition, transmittance is wavelength dependent and, therefore, different for each of the TIRS 
bands. Absorption in the 10.5-12.5 |jm atmospheric window is mainly affected by water vapor that has 
a high spatial variability, since other atmospheric gases, such as CO^ , N2 , and 0^ , are well-mixed, 

and their effect can be considered constant throughout an image for the purpose of this analysis [22,23]. 
The atmospheric water vapor content, at the time of image acquisition, can be obtained from local 
measurements in-situ, or at nearby meteorological stations. The sun photometer measurements of the 
AErosol RObotic NETwork (AERONET), operated by NASA/GSFC [24], are suggested as a good 
source of water vapor data by Qin et al. [18], but they only provide data when a direct line of sight can 
be established between the station and the sun. Thus, sun photometer readings are a good source of 
information only for the processing of daytime and clear sky images. 

Unlike LSE (discussed in the next section), it is not practical to estimate the atmospheric 
transmittance per pixel. Consequently, when using sun photometer point measurements, we make a 
latent assumption that the atmosphere is constant throughout the scene. While this is not always the 
case, this first-order estimation of the atmospheric transmittance for the whole scene, which relies on 
measured water vapor, is preferable to not accounting for it at all [25]. To roughly meet this 
assumption, cloud pixels must be masked. The atmospheric transmittance can then be simulated for the 
entire scene, based on the point measurement of the total content of water vapor in the column and 
standard atmospheric profiles using radiative transfer models, such as Modtran. 

The use of sun photometers for estimating water vapor has some limitations. As mentioned earlier, 
they do not provide information when the sky is cloudy or at night. Furthermore, due to the lack of 
global coverage by ground-based instruments, AERONET data might not be available for all TIRS 
users for estimating the atmospheric water vapor content, especially in the northern latitudes, Asia, 
Australia, and central Africa. In these areas, other means are suggested, such as the use of a total 
column water vapor product provided by European Centre for Medium-Range Weather Forecasts 
(ECMWF) or the Moderate Resolution Imaging Spectroradiometer (MODIS) Precipitable Water 
product (MOD05_L2, MYD05_L2). While these products' spatial resolutions are significantly coarser 
than TIRS, they capture some spatial variability in water vapor distribution. Accordingly, ECMWF 
and MODIS data can be used when the users wish to refrain from assuming constant water vapor 
throughout the scene. The disadvantage of these data products is that in many cases, they do not 
accurately represent the state of the atmosphere at the time of TIRS overpass. Unlike the sun 
photometer that can perform measurements instantaneously or at a close temporal proximity to TIRS, 
the water vapor estimations derived from these products may be produced within several hours of the 
TIRS acquisition, and thus do not capture the exact atmospheric conditions at the acquisition time. 
This is indeed a limitation, and for this reason, these products should be used with caution. As a rule, 
users should select the data source according to its accessibility. When multiple sources of data are 
available, temporal proximity to the TIRS image acquisition and spatial proximity of the water vapor 
measurement site to the study area should be preferred. When facing a dilemma between temporal and 
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spatial precision, the users have to consider this trade-off and make an optimal decision based on 
their experience. 

The results of Modtran 4.0 simulations, conducted for a mid-latitude summer and for a 1976 
standard US atmospheric profile to determine the relation between water vapor and atmospheric 
transmittance, are presented in Table 2. Throughout this paper, the model developed for mid-latitude 
summer is used as an example. Users may opt to use the coefficients for the 1976 standard US 
atmospheric profile where the atmosphere is modeled better. Based on AERONET measurements, the 
water vapor content in the current research area, in the northern Negev Desert, Israel, ranges from 0.5 
to 3 g-cm~^. As can be seen in Table 2, the relation between water vapor and transmittance is close to 
linear. Qin et al. [18] showed that when this relation is evaluated for a larger range of water vapor 
values, it is better to divide the range into several sections and evaluate each of them separately in order to 
achieve a better accuracy. Since the atmosphere in our research area is relatively dry, and 
consequently, the range of water vapor values is relatively small, we are satisfied with the accuracy 
achieved by treating the entire range as a whole. Taking into account a plausible error in the water 
vapor content estimation of 0.2 g-cm , as suggested by Qin et al. (2001) [18] and according to the 
regression coefficients in Table 2, we can determine that the error in the atmospheric transmittance 
estimation is less than 0.031, which is slightly lower than the value obtained by Qin et al. [18]. 

Table 2. Relationship between atmospheric transmittance and water vapor content in the 
column for the water vapor content range of 0.5-3 g-cm . 



Profile 


Estimation Equation 




SEE 


1976US Standard 


rio=-0.1146w-H 1.0286 
Til =-0.1568w-h 1.0083 


0.9882 
0.9947 


0.0094 
0.0086 


Mid-Latitude Summer 


z-io =-0.1134^ + 1.0335 
Til =-0.1546h'+1.0078 


0.986 
0.996 


O.OIOI 
0.0073 



3.2. Determination ofLSE 

The emissivity of land, in contrast to that of the ocean, is significantly different than unity, and 
varies with the heterogeneity of vegetation, surface moisture, roughness, and viewing angle [26]. Since 
LSE can change substantially over short distances, it is important to estimate its value for every pixel 
prior to applying the SWA. Several methods have been suggested to estimate the emissivity for other 
sensors and can also be applied to TIRS. Techniques for emissivity estimations from infrared and 
visible data are reviewed and discussed in detail elsewhere [27,28]. Adapting some of these techniques 
for use with the Landsat-8 requires the use of OLI bands to indirectly estimate the LSE in the TIRS 
bands. For instance, the LSE could be obtained from a land-cover classification, in which the 
emissivity values for each class are assumed [29]. This type of approach is exercised for MODIS LST 
and emissivity products. However, the estimated emissivity in arid and semi-arid areas is potentially 
uncertain, and users are advised to exercise caution in their applications. Of course, since vegetative 
cover tends to change with time [30], good knowledge of the study site and in situ LSE measurements 
of representative ground covers of the different classes that coincide with the satellite overpass are 
desired. However, if the required conditions for the classification approach are not met, it is possible to 
use NDVI by retrieving the proportions of soil and vegetation in order to estimate the LSE [31-33]. 
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4. Sensitivity Analysis 

Several scenarios were considered in order to estimate the possible LST estimation error due to 
misestimating the SWA input parameters: atmospheric transmittance and LSE. These scenarios 
included an LST range of 0-60 °C and a Tio-Tn range from -3 to 3 °C. 

4.1. Sensitivity Analysis to Water Vapor Content 

Since transmittance is derived from atmospheric water vapor content, it is expected that 
transmittance estimation errors will occur simultaneously in TIRS bands 10 and 11. Therefore, the 
sensitivity analysis was conducted for water vapor content, which serves as the input to the model, and 
by which the simultaneously occurring transmittance errors can be estimated according to the 
regression coefficients in Table 2. 

Estimation error of LST is almost independent of temperature change. It changes less than +0.02 °C 
over the temperature range of 0-60 °C, assuming Tjo-Tu = -2.3 (the average case for the Sinai-Negev 
dune field, as seen in Figure 3 A), eio = 0.961, en = 0.971, and underestimating the water vapor content 
in the atmospheric column by 0.2 g-cm for the water vapor content range of 0.5-3 g-cm . This 
minute change is practically neghgible. 

Figure 3. (A) Histogram of Tjg-Tjj computed from the Negev (Israel) — Sinai (Egypt) 
image acquired on 30 June 2013, path 41 row 206; (B) The increase in land surface 
temperature (LST) error (°K) as a result of underestimating the water vapor content in the 
atmospheric column by 0.1 g-cm , for different Tio-Tn scenarios, ew = 0.961 , en = 0.971; 
(C) The increase in LST error (°K) as a result of underestimating the water vapor content 
in the atmospheric column by 0.2 g-cm under the same assumptions as in Figure 3B. 




'^io''^n Water Vapor (gr cm'^) Water Vapor (gr cm'^) 

(A) (B) (C) 

The LST estimation error increases when the atmospheric water vapor content decreases (and, thus, 
the atmospheric transmittance increases). This effect increases rapidly as the brightness temperature 
difference between TIRS bands 10 and 1 1 increases, as seen in Figure 3B,C. However, the contribution 
of the water vapor estimation error to the LST estimation is complex and also depends on the 
emissivity in both channels, as seen in Figure 4. 
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Figure 4. The increase in land surface temperature error (°K) as a result of underestimating 
the water vapor content in the atmospheric column by (A) 0.1 g-cm ; and (B) 0.2 g-cm , 
as a function of the land surface emissivity. Tjo-Tji is kept constant at 1 °K. 
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4.2. Sensitivity Analysis to LSE 
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An error in LSE estimation can occur simultaneously for both of the TIRS bands, but a separate 
error for each of the bands is possible. In the following analysis, we present the example of a 
simultaneous error in both bands. Unlike the relative indifference of the LST estimation error to 
temperature when the water vapor content is misestimated, the LST error is sensitive to temperature 
when the LSE is misestimated. Figure 5A,B presents this sensitivity. 

Figure 5. The increase in land surface temperature (LST) error (°K) as a result of 
underestimating the land surface emissivity (LSE) by (A) 0.001; and (B) 0.005, in both 
bands as a function of LSE and LST. Tjo-Tjj is kept constant at 1 °K, and the water vapor 
content is assumed to be 1.5 g-cm . 
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The LST error increases linearly, with the LST, and decreases linearly when the LSE increases. 
Figure 6A,B presents a similar linear dependence of the LST error by Tio-Tu; as Tio-Tu increases, the 
LST error decreases. In both Figures 5 and 6, the change of the LST error increases at the same rate as 
the LSE estimation error. In addition, the LST estimation error is sensitive to the water vapor content 
when the LSE is misestimated, as depicted in Figure 7. When water vapor content is higher, and the 
atmospheric transmittance is lower, the LST error decreases. Figures 4-7 show that the LST error is 
oppositely related to the LSE. 
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Figure 6. The increase in land surface temperature error (°K) as a result of underestimating 
the land surface emissivity (LSE) by (A) 0.001; and (B) 0.005, in both bands as a function 
of LSE and Tjo-T]]. The water vapor content is assumed to be 1 .5 g-cm and T^o = 3 13. 15 °K. 
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Figure 7. The increase in land surface temperature error (°K) as a result of underestimating 
the land surface emissivity (LSE) by (A) 0.001; and (B) 0.005, in both bands as a function 
of LSE and water vapor content. Tio-Tn is kept constant at 1 °K, and Tw = 313.15 °K. 




0.085 



0.08 



0.075 i. 



0.07 0) 
I- 
co 

0.065 



0.06 



0.92 0.94 0.96 0.98 
Land Surface Emissivity 



(A) 



0.92 



0.94 0.96 
Land Surface Emissivity 



0.98 



(B) 



0.4 _ 




5. Accuracy Assessment of the Proposed SWA 



In order to assess the accuracy of our SWA, we used Modtran 4.0 to simulate the thermal radiance 
reaching the sensor for an input Mid-Latitude Summer atmospheric profile with known LST and LSE. 
The simulated radiance was then converted into brightness temperatures for the TIRS bands, and used 
as inputs to the SWA. The LST estimation errors for 60 different scenarios, featuring different 
combinations of LST, LSE, and atmospheric water vapor content, are presented in Table 3. The RMSE 
for the LST estimation errors in Table 3 is 0.93 °C. 



Table 3. Estimation errors of Land Surface Temperature (LST) for various simulated 
combinations of LST, Land Surface Emissivity (LSE), and atmospheric water vapor content. 



Water Vapor (g-cm ^) 


LST (°C) 


LSE = 0.98 


LSE = 0.97 


LSE = 0.96 


LSE = 0.95 




10 


0.1456 


0.6979 


1.2502 


1.8249 




20 


0.0959 


0.6885 


1.2912 


1.8939 


1 


30 


-0.0021 


0.6410 


1.2840 


1.9372 




40 


-0.0933 


0.5777 


1.2712 


1.9647 




50 


-0.2249 


0.4966 


1.2304 


1.9742 
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Table 3. Cont. 



Water Vapor (g-cm ^) 


LST (°C) 


LSE = 0.98 


LSE = 0.97 


LSE = 0.96 


LSE = 0.95 




10 


-0.3264 


0.2512 


0.7819 


1.3411 




20 


-0.5057 


0.0937 


0.6932 


1.2927 


z 


jO 


—v./a// 


—0.14/9 


0.4919 


1 1 /1 1 o 

1.1418 




40 


-1.1519 


-0.4902 


0.1999 


0.9086 




50 


-1.5968 


-0.8663 


-0.1073 


0.6232 




10 


-1.1184 


-0.6097 


-0.0698 


0.4701 




20 


-0.8608 


-0.2806 


0.3208 


0.9010 


3 


30 


-0.8268 


-0.1852 


0.4565 


1.0770 




40 


-0.9127 


-0.2308 


0.4512 


1.1331 




50 


-1.1507 


-0.4073 


0.3250 


1.0472 



6. Summary 

We presented a SWA for Landsat-8 TIRS data. The ability to derive LST temperatures accurately 
has immediate impacts for potential Landsat-8 users and applications. As data from the new Landsat-8 
are distributed freely, its high resolution thermal abilities can allow for new scientific advances in earth 
science. Although the spatial resolution of TIRS is degraded to 100 m, in comparison with the 60 m 
resolution of the Landsat-7 ETM+ thermal band, the added value of having two bands is more accurate 
LST estimations with TIRS than with its predecessor. In addition, the 100 m resolution is sufficient for 
water consumption measurements over fields irrigated by center pivot systems [9], as well as for other 
uses over relatively homogeneous areas. 

In this paper, we have outlined possible strategies for evaluating the atmospheric water vapor 
content and LSE. When selecting a data source for estimating the water vapor content, the tradeoff 
between temporal and spatial precision has to be considered by the users on a case-by-case basis. 
Strategies to estimate LSE, based on techniques employed by previous sensors, have been suggested. 
These techniques can still be refined and adapted for OLI's new spectral band configuration, which is 
different than its predecessors. Therefore, future work should focus on evaluating the best methods of 
estimating the SWA input parameters (e.g., atmospheric transmittance, LSE), the cross-sensor 
validation of the input parameters and the resulting LST. 
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